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Abstract 

In this paper we present the calculation of a scalar pentagon integral with 
two consecutive massive external legs having an equal mass propagator 
embedded between them. We also deal with the two situations where the 
farest external leg is either massive or not. The relevance of the calculation 
comes from its application in many perturbative QCD calculations as well 
as in QCD corrections for weak precesses. 
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I. INTRODUCTION 


A series of rare elementary processes involving more than two particles in the final 
state are going to be measured with increasing precision. The multiplicity of the final 
state makes it difficult to extract predictions by the standard gauge theories even if 
semplihcations arise when either participants are all massless or only some of the external 
particles are massive. However more accurate rate measurements of processes with heavy 
quark hadrons in the final state will soon be available as is the case of the CHORUS 
experiment where direct evidence for the associate charm production in charged current 
neutrino nucleon scattering has been shown [1], In the one loop calculation of such 
processes we encounter pentagon integrals with a massive line as skeched in figure 1 where 
massive particles are bold, massless ones thin and dashed ones can be either massive or 
not. 



In general the inclusion of masses makes things more involved, although the calcula¬ 
tion simplifies when either external masses are equal to each other or they are equal to the 
internal masses or both eventualities occur as it is often the case in normal gauge theories. 
Recently a lot of progress has been made in the technics for perturbative calculations with 
different approaches. A non-comprehensive list is given in [2] and reference therein and 
in [3], [4], [5], [6] and [7]. In particular adopting the dimensional regularization approach 
for Feynman parametrized integrand the authors of ref. [3] derived simplifications and 
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recursion formulas by the implementation of algebraic technic. Using these methods the 
problem of the evaluation of a one loop n points scalar integral is translated to the eval¬ 
uation of a combination of n — 1 points scalar integrals and the original n points integral 
in D = 6 — 2 e dimensions; moreover the original n points one loop integral can be repre¬ 
sented as the solution of a partial differential equation system. In the present paper we 
use this approach to perform the calculation of the pentagon integral represented in figure 
1. Other massive pentagon integrals have been recently evaluated in next to leading order 
calculations of processes in which an Higgs particle can be generated at hadron collid¬ 
ers. In particular two independent groups report the NLO corrections for the process in 
which an Higgs particle is generated together with a tt pair, [8] and [9]. Another NLO 
calculation involving massive pentagon integrals is given in [10] in which the final state 
considered consists of an Higgs particle plus two jets. The general methods employed here 
do not concern with the specific processes and the results must be analitically continued 
to describe a specific process. Finally only the most simple tensor integral is given while 
we postpone other cases to a dedicated paper [11], The paper is organized as follow: in 
section II relevant formulas from ref. [3] are collected, in section III they are applied to the 
scalar massive pentagon represented in figure 1 transforming it in a combination of four 
points integrals; section IV is devoted to four point integrals evaluation and in section V 
more simple tensor integral (vector) is given with the conclusions. The initial condition 
for the differential equations originated in the four points evaluation are calculated in the 
appendix. 


II. BASIC FORMULAS 


The starting point is the integral in D — 4 — 2e dimensions 
d 4 ~ 2£ l 1 


Tj 1 /i 


2s 


( 27 t ) 4_2£ (/ 2 - M 2 )((/ - Pl y - Mf )...((/ - p n _ 0 2 - Ml) 
with the momenta hi taken to be outgoing, k 2 = mf and 


( 1 ) 


Pi = k i 

3 = 1 
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Po = Pn = 0 . 


( 2 ) 


Applying Feynman parametrization, Wick rotating and integrating over loop momentum 
this integral can be cast in the form 


= i (—l) n (4vr) £_2 fi 2e I> 


I n = T(n — 2 + e) [ d n ai 5( 1 - V. 

Jo 


D(di 


in-2+e 


with 


= '£,s, 

i,j =1 


iy d i Clj 


and the matrix S given by 


with pa = 0 and 




( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


Pij = ki + k l+1 + ... + kj- 1 (7) 

for i < j. We will not repeat the derivations obtained in ref. [3] but, to introduce notation 
and to be self-consistent, in the rest of this section we just collect relevant formulas 
that will be used in section Ill and IV. Performing a projective transformation [12] with 
parameters a t in such a way that the denominator in Eq.(4) has no dependence the 
definition of a new matrix follows (indices are not summed) 

Pij 'Jiy ^i A?* (3) 


Using the following definitions 



N n = 2 n ~ 1 det(p ) 
n = N n p~ l 


7 i — } " Pij a j 

3 =1 


A n ^ ^ Tjij Oij Oij 
i,j =1 


( 9 ) 
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the authors of ref. [3] find 


In 

1 di n 

n — 4 + 2e dcti 


1 

2N~ n 

1 

2iVn 


^ li In-i + {n — 5 + 2e) A n J,f 6 2e 
_ i=l 
n 

^ ^-1 + ( n - 5 + 2e ) 7* in = ^' 2£ 

-j= 1 


( 10 ) 

(li) 


where { stands for the n — 1 integral with the denominator obtained from an I n integral 
eliminating the propagator between legs 1 — 1 and i; once Feynman parameter has been 
introduced in the usual way for I n the denominator in , is obtained putting a* = 0. By 
the observation that I 4 and 1$ are finite in 6 dimensions, performing one-loop calculation 
one can limit to evaluate 


h 

dl 4 

dai 



( 12 ) 

(13) 


taking only the divergent part from the L 


U) 


integrals in Eq.(13). 


III. PENTAGON WITH ONE MASSIVE LINE 


To write down the integral in figure 1 we set k\ — k\ — Mf = Mf = Mj = M| = 0, 
kl = kl = Mi = m 2 and = q 2 giving 

d 4 ~ 2 H 1 


X5 /1 


2e 


(2tt) 4 - 2£ ( l 2 )(l - Pi) 2 ((l - P 2) 2 - m 2 )(l - p 3 ) 2 {l - p A ) 2 


(14) 


and 


h = r(3 + e) [ d 5 ai 5( 1 -V, 

J 0 


V(a. 


.^3+6: 


(15) 


with T> given in Eq.(5) and the matrix S given by 
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s 


( 0 

0 

m 2 - S 12 
—< S 45 

V -q 2 


0 

0 

0 

— ^23 
— S51 


m 2 - S 12 


0 

2 m 2 
0 

m 2 - s 34 


-S 45 -q 2 \ 

—S23 — S51 

0 m 2 — S 34 

0 0 

0 0 / 


(16) 


with = (ki + ki + \) 2 . We define = Sj,i+i — m 2 and in the following we will 

assume S12, S23, S34, S45, S51, q 2 < 0. Performing the projective transformation with 


Ot\ — 

0,2 — 

OL 3 = 

0^4 — 

CD 5 = 

we get for the Pij matrix in Eq.( 8 ) 


-S23 s 34 


s 45 S 51 s 12 

S45 ^34 

-S23 S 51 S12 
<S45 S 51 
^23 ^34 S12 
$51 Sl2 
-S 23 S34 S45 
Si2 S23 
S34 S45 S51 


P= 2 


/0 0 1 1 

0 0 0 1 

1 0 2 M 2 0 

1 1 0 

\ A 1 1 


with 


A = 


q 2 S 23 
S 45 S 51 


A\ 
1 
1 


0 0 
00 / 


(17) 


(18) 


(19) 


M 2 = - 


m S 45 S 51 
S12 S23 S34 


( 20 ) 


If k\ = f = 0 we only have to take A = 0 in Eq.(18). The coefficient relevant for the 
evaluation of the pentagon by Eq.(12) are given in the table 1, keeping apart the case 
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q 2 = 0. Due to the presence of masses we have not cyclic relations between the coefficients 
but only the relations 


74 = 721 


Cx. 4<—>C*!2j Q:5<—>Q!1 


75 = 71 


OL 4<—»Oi2, 


( 21 ) 


Table 1. Coefficients to be used in Eq.(12) 


par 

any q 2 

to 

II 

O 

n 5 

1 + M 2 - (l + 2 Af 2 ) A + Af 2 A 2 

1 + Af 2 

7 i 

«1 — OL 2 + O3 — A «3 + O4 + 2 Af 2 04 

— 2 Af 2 A 04 — 05 — 2 Af 2 05 + 2 Af^ A 05 

OLI — OL 2 H - £^3 + (l + 2 M 2 ) (074 — 0^5) 

72 

—+ 0^2 — + A 0^3 + 0^4 — 2 A OL 4 — 2 M 2 A 0^4 

+2 M 2 A 2 o;4 + a?5 + 2 M 2 a5 — 2 M 2 A 0^5 

—01 + 02 — 03 + 04 + (1 + 2 Af 2 )o5 

73 

cki — A c^i — Q2 + A 0L2 — 2 A ol 3 

+A 2 OL 2 > — C ^4 A 0^4 + CK5 — A 0^5 

Ol — 02 + 03 — 04 + 05 

74 

QL\ + 2 M 2 a\ — 2 M 2 A o^i + a2 — 2 A «2 

—2 .A/f 2 A c^2 H - 2 M 2 A 2 oi2 — 0^3 H - A 0^3 -K 074 — 0^5 

(1 + 2 Af 2 ) Ol + 02 — 03 + 04 — 05 

75 

— ai — 2 Af 2 ai + 2 Af 2 A ai + 02 + 2 Af 2 02 

— 2 Af 2 A 02 + O3 — A O3 — 04 + 05 

(1 + 2 Af 2 ) ( a 2 - 01) + 03-04 + 05 


In terms of new kinematical variables cti, A and M the denominator in the / 5 integral 
represented in Eq.(5) is given by 


M 2 d3 2 0\ O3 Q’i O4 CI2 CL4 A Oi O5 0-2 CL5 O3 O5 

CL 3^ CL\ OI3 CL\ OL\ OI2 Ol 4 CK4 OL 5 Q?2 CL 5 CK3 O5 


( 22 ) 


and the four points denominators in the ] integrals in Eq.(12) can be obtained putting 
«j to zero in the expression above. It is easy to verify the relations 
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( 23 ) 


/,[ ) (Q!i, «2j « 3 > « 5 ) 
a 2, «3> <* 4 ) 


I4 ^(cKi, «3, 0(4, CX5) 
l[ \&2, OI3, CK4, CK5) 


0:4—>0:2, Q!5-<-^Q;i 


CK5— >& 1 , Ot 4 <—>0:2 


In the next section we proceed to the evaluation of \ /j 2 ^ and lf } using the set of 
partial differential Eqs.(13). 


IV. FOUR POINTS FUNCTIONS 

Here we evaluate the integrals and \ corrsponding to massive boxes with an 
internal massive line, in the variables defined in Eqs.(17, 19, and 20) and translate the 
integrals I 4 that are well known and correspond to massive boxes with massles internal 
lines. 


A. The integral / 


(i) 


After putting a\ = 0 in Eq.(22) we have for the denominator in 

Iff 2 03^ 0,2 O4 O2 0,5 CL 3 O5 

Ol 3^ CK2 OL4 CC 2 C»5 CK3 CK5 

Before solving the integral we perform the following kinematic transformation: 


(24) 


0.2 = M c 2 
ck 3 = M c 3 
ol 4 = C4/M 
«5 = C 5 /M. 


(25) 


In terms of the new variables we get for the denominator: 

O 3^ (Z 2 O3 

C3 2 C 2 C4 C 2 C5 C3 C5 


( 26 ) 


and 



N 4 = 1/2 

/O 0 


0 0-11 


1 - 1-2 2 


.0 1 2 - 2 , 


The only divergent three points functions extracted by the expression above are A, and 
I ?j obtained putting 02 = 0 and 03 = 0 respectively; these correspond to two two-mass 
triangles, while the other two obtained putting a 4 = 0 and a 5 = 0 respectively are three- 
mass triangles checked to be finite. At the Oiy£~ l ) we have 


/W = r(l+e) 


1 log( c s) 

2 £ 2 C 4 £ C 4 


+ O{e 0 ) 


f(2) _ r(l + e) , (C5 

3 _ £(C4-C 5 ) g U 


+ 0(£°). 


The system of partial differential equations in Eq.(13) is then given by 


dlj’ T(l + e) [1 


C 4 - C 5 \£ 


- 1 - — + 2 log(c 4 ) -2 — log(c 5 ) 


dl^ 2 T(1 + e) 


C4 C5 


log — 


with the solution 


I 4 = r(l +£) - log(c 4 ) + log ( 1 - — ) - log(c 4 - c 5 ) ) +2 log(c 4 ) log(c 5 - c 4 ) + 


2 log(c 5 ) log(c 4 — C5) + k\ . 


The integration constant k\ is evaluated in the appendix and its value is 


fci = + 5C(2) 


where LE is the dilogarithm function and £(2) = 7 t 2 /6. Afetr some manipulation we have 
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r ( i + £ )|g-v(| 


2 Li 2 ( 1 — — j + y- 
c 5 / 6 


Reintroducing the original variables inverting Eq.(25) we get for I\ 


(i). 


/| 1} = r(l + e) 


[ _ log * fsa) - 


OL 4 


2 Li 2 1-— + — 


71 


y 2e 2 ‘ ° \a 5 J V* «5 ) ’ 6 

Being l\ ]) independent from A its value does not change in the limit q 2 —> 0. 


(33) 


(34) 


B. The integral 

Here and in the following subsection we proceed performing the same steps as in the 
derivation of The limit q 2 —> 0 now gives a different situation; in fact in this limit 
there will be three divergent three-point integrals extracted by I\ ; so as explained in [3] 
the limit procedure is not smooth and the two case have to be taken separately. 


1. I 


( 2 ) 


q 2 + o 


In this case I\ is a three external mass box but, differently from 
external masses different from each other and so it needs evaluation. After 
in Eq.(22) we have for the denominator in I\ J 

M 2 a , 3 2 cli a 3 a, 1 04 A cq 05 03 05 

OL 3 2 OL 1 OL 3 Oil OI4 OL\ O5 O3 O5 


it has all 

putting a -2 = 0 

(35) 


Rescaling the variables with 


= C\/M 
a 3 = Mc 3 
a 4 = M c 4 
a 5 = c 5 /M 

A = 8 /M 2 (36) 

we get 
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(37) 


giving JV 4 = 1/2 and 


a-3 2 «i 03 Ol 5 «1 «5 «5 

C3 2 Cl C 3 Cl C 4 Cl c 5 c 3 c 5 


0 0 


0 0 


-5 25(1 - 5) 25-1 
1 25-1 -2 . 


The only divergent three points functions extracted by the expression above are il and 
I' s obtained putting op = 0 and 03 = 0 respectively; these correspond to two two-mass 
triangle, while the other two obtained putting a 4 = 0 and <25 = 0 respectively are three- 
mass triangle checked to be finite. At the (9(e _1 ) we have 

f( l ) _ TVi 1 1 1 ^ 1 /orO 


I3 = r(l+e) 


2 £ 2 C 4 


+ 0(e u ) 




The system in Eq.(13) is then given by 


= £(i+£) ri ( c jl _ s \ + 2 c Jllog(cs) _ 25log(5C4) 

OC4 C5 — 5 C4 [£ \C 4 ) C4 

d¥p = 2 T (1 + e) /5 q\ 
dc 5 c 5 - 5 c 4 V c 5 / 

with the solution 

iP = T(1 + e) y log(c 4 ) + 2 log(c 4 ) log(c 5 ) - 2 log(c 5 - 5 c 4 ) log ^1 - 
+2 Li 2 + lo g 2 ( c 5 - 5 c 4 ) - 2 log(c 5 ) log(c 5 - 5 c 4 ) 

+ log 2 f 1 - + 2 log(5) log(c 5 - 5 c 4 ) + 2 log(c 4 ) log f 1 - + k 2 
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The integration constant fc 2 is evaluated in the appendix and its value is 


h = 2^2 + 7 lo g( 5 ) _ \ lo S 2 ( 5 ) “ Li 2 (l - 5) - 2 C(2). 


(44) 


After some maipulation we have 


if = r(i + e) 


£ 2 £ J 2 2 V c 5 / V C 5 


-/o') 

Reintroducing the original variables inverting Eq.(36) we get for I\ : 


/f = T (1 + e) 


1 1 
“77 -hi 


£ 

2 


(M A a 4 ) 


2 e 


1 


Li 2 (l — M A) — 2 Li 2 1 


2 

A 04 

OL 5 


- -log (M A) - log 

2 \ a 5 


A ol 4 


(45) 


(46) 

(47) 


S. if, g 2 — 0 


^(2) 

In this case I\ is a two external mass box. Putting A = 0 the denominator in 
Eq.(35) became 

M 2 a 3 2 ai d 3 ai 04 03 d$ 


+ 


+ 


+ 


(48) 


cx 3 ot 1CO3 ot 1 04 CK3 05 

Rescaling the variables as in Eq.(36) we get 

2 

tt.3 Ol 03 Ol O4 ^ 03 05 
c 3 2 C 1 c 3 C 1 c 4 C 3 C 5 

The divergent three point functions if, if and if are obtained putting d\, a 2 and 
d 3 = 0 respectively. 


if = r(l + e) ( -b + -L log(c 5 ) ) + O(e 0 ) 


if = r(i + £) 
if = r(i + 0 


C4 £ C4 

~pnr + tt - lo s( Ci C4 )) + c ’( £ °) 

S C5 ^ C5 


e(ci-c 5 ) V c i 


log (-) + 0 ( 0 ) 


(50) 

(51) 

(52) 
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while is a three-mass triangle checked to be finite. The partial differential equation 
system is given by 

d /| 2) 


k = affi±£> log hd 

c 1 Cl — C5 V C 1/ 


dc- 


di? ] 


dc?, 


= 0 


dl 


( 2 ) 


r (i + «) (1 


Dca 


c 4 


+ 2 log(c 5 ) 


dl 


( 2 ) 


dc§ C 5 

with the solution 

/| 2) = r(l + e) 


2 T(1 + e) /1 Ci , . 

' +-log(ci) - 


C 5 


£ Ci - C 5 


Cl — c 5 


■ log(c 5 ) + log(c 4 ) 


(53) 


Cl 


(log(c 4 ) + 2 log(c 5 )) - log 1 -- log(ci - c 5 ) - 


C 5 


Cl 


2 log(ci) log 1 - +2 log(c 4 Ci - C4C5) log(c 5 ) 

V C5 1 

Cl' 


—2Li 2 J + k 

The integration constant k :i is evaluated in the appendix and its value is 

3 


*3 = ^5 - 8C(2) 


After some manipulation we find 
/f =T( 1 + e 


+ - log 


2e 2 




(— V 

V c 5/ 


Reintroducing the original variables we have 

\ 2 e 


/| 2) = T (1 + e) 


04 ^ 2 £ (Mq 5 ) 


M) 2e 2 




log 


+ 2 Li 2 


CK 4 

M 2 a 5 


Cl 

c 5 


+ 2Li 2 


4C(2) 


CH 5 


(54) 


(55) 


(56) 


4 C (2) 
(57) 


C. The integral I ® 

Putting 03 = 0 in Eq.(22) we eliminate the massive propagator and obtain the easy 
(opposite) two mass box [13] or the one external massive box if we take respectively q 2 7 ^ 0 
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or q 2 = 0. These integrals are well-known and are reported also in [3]. Here we just put 
these integrals in the kinematics specified in section 3. 


1. if, q 2 ± 0 


After putting 03 = 0 in Eq.(22) the denominator is given by 

(X 1 CLa do CL,a A Oi Or O9 Or . . 

-1-f-- H-- (58) 

0\ 0^4 0^2 O4 0\ 0^5 O2 O5 

Using Eq.(4.44) from the third paper in ref. [3] the integral in the kinematics of section 3 
reads 

(3, _ 2T(1 +£) 


TW _ 
J 4 — 


1 -A 


/ (ai o 4 1' - (a 2 a 4 f + (a 2 o-,r - ( — j * 5 


- - log 2 ( 

2 y cx\ ex 4 


—Li 2 ( 1 — —4 — Li 2 fl-—") + Li 2 (l — A) — Li 2 (l-^) - Li 2 (4 - — 
a 2 J V a 4 J \ o.i ) \ a 5 


2. if, q 2 = 0 


Putting A = 0 in Eq.(58) the denominator of this integrals is given by 

Oi O4 o 2 O4 o 2 O5 


(60) 

Ol\ CY4 CK 2 CO4 0 L 2 OL 5 

Using Eqs.(4.27, 4.40) from the third paper in ref. [3] the integral in the kinematics of 
section 3 reads 


if = r(l + e) 


4 ((oq a 4 ) £ - (a 2 a 4 y + (a 2 a b ) £ ) - log 2 ( 

a 2 a 5 


—2 Li 2 ( 1 — — j — 2 Li 2 (1 — — ) -4C(2) 
a 2 J \ a 4 1 


(61) 


V. CONCLUSIONS 

An expression for the scalar pentagon integral shown in figure 1 can be built via 
Eqs.(12), (23), the four points integrals evaluated in the last section and the coefficients 
in table 1. The expressions for are very long and are not reported. 


(59) 
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More familiar kinematics is realized by replacing the variables cq, A and M 2 with 
their definitions in terms of Sij , q 2 and m 2 . Tensor integrals will be considered in a separate 
paper [11], however the simplest one of them, the vector integral, is related to the scalar 
integrals with one Feynman parameter in the numerator by the following relation 

4 D n -* Fm (62) 

in which the arrow means integration over loop momentum /, the integrand numerator is 
in the square brackets and 

n— 1 

'P 1 ' = a i+'Pi ( 63 ) 

i= 1 

with pi given in Eqs.(2). The integrals h[dj\ can be evaluated by [3,14] 

= + o(e) ( 64 ) 

2iVs tr 

where p dehned in Eq.(9) is deduced by p given in Eq.(18) 


/ 

1 


-1 

1 - A 1 - 2M 2 (-1 + A) 

-1 + 2 M 2 (-1 + A)\ 


-1 


1 

— 1 + A 1 + 2 (—1 + M 2 (—1 + A)) A 

1 — 2M 2 (—1 + A) 


1 - A 


— 1 + A 

(-1 + A) 2 —1 + A 

1 — A 


1 - 2 M 2 (- 

1 + A) 

1 + 2 (-1 + M 2 (-1 + A)) A 

-1 + A 1 

-1 

V 

-1 + 2M 2 (- 

-1 + A) 

1 — 2M 2 (—1 + A) 

1 - A -1 

1 ) 


Higher tensor integrals can be evaluated considering that they are linked to scalar integrals 
with more powers of Feynman parameters in the numerator [3]. Such a decompositioncan 
can also be organized in a way that drastically reduces numerical instabilities genarated 
by the presence of inverse powers of Gram determinants [14]. Besides the deep inelastic 
case mentioned in the introduction, the results obtained in the present paper with q 2 ^ 0 
can be useful in the evaluation at one loop of the decay amplitude of a real W boson or 
a virtual photon in a heavy quark-antiquark pair and two light quarks. Let us consider 
the case in which all massless particles and k§ are gluons, then the pentagon studied with 
q 2 = 0 can be identified with one of the four pentagon in the perturbative evaluation 
of the one-loop associated production of heavy quark in the gluon-gluon-fusion with a 
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gluon in the final state (gg —> QQg ); in this case, indeed, pentagons are found in which 
the propagators form chains with 1, 2, 3 and 4 equal mass fermions the first of which is 
calculated in the present paper while the other ones can be calculated analogously. 
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APPENDIX A: FOUR POINTS INITIAL CONDITIONS 

In this appendix we report the calculation of the integration constants for the four 
points integral of section four systematically neglecting O(e) terms. Instead of report¬ 
ing all length passages, we give the steps that can be followed by programs of function 
manipulation like the used Mathematica. 


1. Integration constant for ^ integral 


The point chosen to evaluate is C 2 = C 3 = ( 04 / 2 ) = C 5 = 1 where the expression 
in Eq.(31) gives 


i? = r(l + e) ( j log( 2 ) — 3 C( 2 ) + Aq 


(Al) 


The expression for the integral at the point selected deduced using Eq. (4) and the first of 
Eqs.(9) is 


(i) _ r (2 +£) .5(1-E,» 


J 4 — 


' drai —;— . . 

0 (O3 + 2 ® 2®4 T U 2 a 5 T O3O5) 


(A2) 


The factor 2 is given by IIq. Renaming 02 with x, <23 with y and <24 with z, and performing 
the transformation x —>1 — x, y —> x — y and z ^ z we arrive at the expression 
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/(.) = ip ±£) r dx r dy r dz 

2 Jo Jo Jo 


(x 2 — 2 xy — z + y z + y — \ z) 2+£ 


Putting apart the Gamma function for the moment, the z integration gives 


— (x 2 + y — 2xy) 1 £ 
(1 + e) (l + x-2y) 


2 1+£ (2x 2 + y-hxy + 2y 2 ) 
(1 + e) (l + x-2y) 


2 \ — 1 — £ 


The two integrals can be evaluated by shifting both in y 


y~y + -_ 


simplifying the x integral 


2 1+£ (x + 2y - Axy) 1 £ 

(1 + e) (1-2 y) 


2 2+2e {x + 2 y — 6 xy + 4 y 2 ) 
(1 + e) (1-2 y) 


2\-l-e 


and inverting the integration order [12]. After some manipulation and expanding some 
Hypergeometric and Generalized Hypergeometric functions the result is 

W = r (2 + <r) (Tj - T + I log( 2 ) + i + 2 <( 2 ) - log( 2 )) . (A 7 ) 

Finally, making the substitution T(2 + e) — (1 + e) T(1 + e) in Eq.(A7) and taking into 
account Eq.(Al) we End k\ in Eq.(32). 


2. Integration constant for integral with q 2 ^ 0 


To evaluate the integration constant we evaluate the integral in the point: 


5-1 

1 

5-1 

1 

1-5 


c 4 = —1 
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lie, = 


(A10) 


(1 -<5) 3 

The integrations are trivial but the expression is very long. The q chosen cannot be 
simultaneously positive so we checked the result in the point (5/2) = ci = c 3 = c 4 = c 5 = 1 
where the expression in Eq.(43) gives 

/f =T( 2 + e) 


_L _ i + I log( 2 ) + | + §C( 2 )-l°g( 2 ) - \ log 2 ( 2 ) 


(All) 


The expression for the integral at the point selected is 


/| 2) = T(2 + e) 


d 4 a,i- 


5(! - E,; 


(A12) 


o (03 T Q4O3 T u 1 04 T 2 cq(Z 5 T (2305)^^ 

Renaming a, with x, 03 with y and 04 with z, and performing the transformation x — > 1— x, 
y — > x — y and z ^ z we arrive at the expression 

if' 1 = T(2 + e) / dx I dy I dz- 


'■i r x rv X 


jo jo Jo (-2 xy+ yz + x + y - z) 2+£ 

Putting apart the Gamma function for the moment, the z integration gives 


(A13) 


(x + y — 2xy) 


— 1 —£ 


+ 


(x — 2xy + y 2 ) 


2 \—1— £ 


(A14) 


(1 +£) (1 -y) (1 + e) (1 - y) 

Performing the x integration before and adding and subtracting terms we find Eq.(All). 


3. Integration constant for if' 1 integral with q 2 = 0 

The point chosen to evaluate I\ is given by 2ci = (C 3 / 2 ) = (C 4 / 2 ) = C 5 = 1 in 
which the expression in Eq.(54) gives 


/f = r(l + e) log( 2 ) - log 2 ( 2 ) + 5 C(2) + k 3 . 
The expression for the integral at the point selected is 


(A15) 


(2) _r (2 + £) f l M 


rv_ 
J 4 — 


d 4 a,i 


(A16) 


2 Jo + O1O3 + O1O4 + ^ 0,3(1 f 

Renaming a 3 with x, 04 with y and 05 with z, and performing the transformation x —> 1— x, 
y — > x — y and z —> z we arrive at the expression 
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?( 2 ) _ 
J 4 — 


r(2 + s) 


dx dy dz 


(\x 2 - y 2 - \xz + yz - \x + y - \z + |) 


i\ 2 + £ 


. (A17) 


Putting apart the Gamma function for the moment, the z integration gives 
2 2+2e (1 - x) _1 “ £ (1 - x + 2 y)~ l ~ £ 2 2+2e (l-2x + x 2 + Ay-A y 2 )~ l ~ £ 


(1 +e) (1 + x - 2 y) 
Performing the shift 


(1 + e) (1 +x - 2y) 


(A18) 


y^y + 


X 


(A19) 


gives 

2 2+2e (1 - x)- l ~ £ (1 + 2 y)- l ~ £ 2 2+2e (1 + Ay-Axy -4 y 2 )^ 

(1 + 5 ) (1-2 y) (1 + 5) (1 — 2 y) ’ 1 J 

finally performing the x integration before and adding and subtracting terms before ex¬ 
panding in 5 the result is 

/f = r(2 + e) (+j - + + 1 log( 2 ) + | - 3 C(2) - log(2) - log 2 ( 2 )) . (A 21 ) 

Substituting T(2 + e) = (1 + e) T(1 + e) in Eq.(A21) and taking into account Eq.(Al5) 
we find in Eq.(55). 
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